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We construct the general spherically symmetric and self-similar solution of the Einstein-Vlasov 
system (collisionless matter coupled to general relativity) with massless particles, under certain 
regularity conditions. Such solutions have a curvature singularity by construction, and their initial 
data on a Cauchy surface to the past of the singularity can be chosen to have compact support in 
momentum space. They can also be truncated at large radius so that they have compact support in 
space, while retaining self-similarity in a central region that includes the singularity. However, the 
Vlasov distribution function can not be bounded. As a simpler illustration of our techniques and 
notation we also construct the general spherically symmetric and static solution, for both massive 
and massless particles. 


I. INTRODUCTION 


In an astrophysical process of gravitational collapse of 
matter there is always some physical effect opposing the 
force of gravity, so that not all the matter necessarily 
ends up in the final black hole. Since the initial work of 
Choptuik [jij, it has become clear that in many physical 
systems it is possible to fine tune the initial conditions 
so that collapse can result in black holes of arbitrarily 
small mass. This gives rise to interesting effects which 
are now studied as a particular branch (called type II) 
of Critical Phenomena in Gravitational Collapse (see || 
for a review). In a perfect fluid it is the pressure that 
opposes collapse. In field theories, it appears to be the 
fact that the matter field propagates at the speed of light. 

Here we study another type of matter, a cloud of point 
particles that do not interact directly through collisions 
but only through their averaged gravitational field. This 
system is described by the coupled Einstein and Vlasov 
equations. In collisionless matter it is the tendency of 
particles in the cloud to miss each other that acts as the 
source of dispersion. In a spherically symmetric situa¬ 
tion, such as we shall study here, one can also think of 
their angular momentum as the force opposing gravity. 
(In spherical symmetry, each particle has angular mo¬ 
mentum, but the total angular momentum of all particles 
is zero.) 

In numerical simulations of the gravitational collapse 
of asymptotically flat, spherically symmetric collision¬ 
less matter configurations Rein, Rendall and Schaeffer 
|jj have not found any sign of critical phenomena at the 
black hole threshold. Olabarrieta and Choptuik ||] 
have not found type II critical phenomena either, but 
have found some evidence of type I critical phenomena, 
where a metastable static solution acts as an interme¬ 
diate attractor at the black hole threshold. They have 
found the expected scaling of the lifetime of the interme¬ 
diate attractor with distance to the black hole threshold. 
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The critical exponent, however, varies from 5.0 to 5.9 be¬ 
tween families of initial data, for a quoted uncertainty 
of 0.2. The Vlasov distribution function of the critical 
solution was not found to be universal, and the metric of 
the critical solution was found to be universal only up to 
a rescaling of space. 

We want to explore the problem of existence of critical 
phenomena in the Einstein-Vlasov system by looking for 
both static and continuously self-similar (from now on 
CSS) solutions of the system and then testing for their 
stability against linear perturbations. A static solution 
with exactly one growing mode is a potential critical so¬ 
lution for type I critical collapse, and a self-similar solu¬ 
tion with exactly one growing mode is a potential critical 
solution for type II critical collapse. In this paper, we 
construct the general static and the general self-similar 
solutions, assuming some regularity conditions. In order 
to avoid introducing a scale, we further restrict ourselves 
to massless particles in the self-similar case. Their per¬ 
turbations and the extension to massive particles will be 
studied elsewhere. The static solutions of the Einstein- 
Vlasov system have been extensively studied already. We 
deal with them here partly as an introduction to the 
methods we use later in the CSS case but also for the 
light they shed on type I critical phenomena. To our 
knowledge little is known about CSS solutions: The CSS 
equations were derived in |S), but no CSS solutions were 
found. Spherically symmetric CSS solutions for rotat¬ 
ing dust, where L(r) can be specified freely in the initial 
data on top of p(r) and v(r), were constructed recently 
by Gair @]. 

We begin by fixing notation and stating the field equa¬ 
tions for the Einstein-Vlasov system in spherical symme¬ 
try. We then introduce the idea that the general solution 
of the Vlasov equation on a fixed spacetime can be given 
in terms of an arbitrary function of three variables that 
are conserved along trajectories and which can be given 
as quadratures if the spacetime has an additional con¬ 
tinuous symmetry, besides spherical symmetry. This is 
true in particular for both massive and massless particles 
in static spacetime, and for massless particles in a CSS 
spacetime. 

The further requirement that the resulting stress- 
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energy of the collisionless matter is consistent with the 
assumption of staticity or CSS reduces the free function 
of three variables to one of two variables. We show this 
first for the static case, and discuss an interesting sub¬ 
tlety. We also construct some explicit solutions numeri¬ 
cally. This section can be omitted by the reader who is 
only interested in the self-similar case, but we feel that 
it serves as a useful warmup. 

In discussing the CSS case, we first discuss the con¬ 
served quantities, and classify the possible particle tra¬ 
jectories on a CSS spacetime. We then find the general 
solution on a fixed background spacetime, and reduce it 
to the general self-gravitating solution. We again con¬ 
struct some explicit solutions numerically, and discuss 
the general features of these solutions. The implications 
of our results for critical phenomena in gravitational col¬ 
lapse are discussed in the Conclusions section. 


II. THE EINSTEIN-VLASOV SYSTEM IN 
SPHERICAL SYMMETRY 

A. Collisionless matter 

In this section we briefly review the essential ideas 
and equations needed in the rest of the Article. For a 
longer introduction to the Einstein-Vlasov system see || . 
We add some considerations about the integration of the 
Vlasov equation and its relation to the symmetries of the 
spacetime. 

The Einstein-Vlasov system describes the evolution of 
a statistical ensemble of non-interacting particles cou¬ 
pled to gravity through their averaged stress-energy. It is 
possible to work with particles of different masses, but in 
this work we restrict ourselves to a single particle species, 
with mass to > 0. We keep to arbitrary, so that all the 
equations in this section are valid both for massless and 
massive particles. 

As in classical statistical mechanics we describe the 
state of a many body system with a positive distribution 
function over the phase space of the system. As the parti¬ 
cles do not interact except through their averaged stress- 
energy, we can use a distribution function /(aT, p v ) on 
the one-particle phase space, where x^ are coordinates 
in spacetime and the coordinates jp are constructed us¬ 
ing the metric at the point x^ from the coordinates 
in its cotangent space at that point. As we always have 
guvP^p 1 ' = —to 2 and p° > 0, we can consider the distri¬ 
bution / as a function f(x^,p l ). Greek indices denote 
the range 0-3 and Latin indices the range 1-3. 

The free-fall trajectories of the particles define a con¬ 
gruence of curves on phase space, tangent to the Liouville 
operator (or “geodesic spray”) 


time s through ds = mda for massive particles. This 
allows us to discuss both massive and massless particles 
simultaneously. As the particles do not interact directly, 
the evolution of / is governed by the Vlasov equation 

Cf = 0. (2) 

The stress-energy tensor at each point is obtained 

by integrating the distribution function over momentum 
space: 

T^(x) = f —yf z gf{x 1 p)p lJj p 1 ', (3) 

JP(x) ~P0 

where P(x) is the 3-momentum space at the point x^ 
and po is determined from p 1 and the metric. The Vlasov 
equation is a sufficient condition for stress-energy conser¬ 
vation. The particle number current given by 

N^(x) = f —V^gf{x,p)p 11 . ( 4 ) 

Jp{x) ~P o 

is also conserved. Note that the distribution function 
f(x,p) is a scalar on phase space even though there is a 
factor \/—g in the integrals. The natural volume mea¬ 
sure on the phase space is — dx 0123 A dpom, which can 
be rewritten as dx 0123 ) A dp 0123 ). It is the 

second factor that appears in the integrals (|J) and 
([|), which are themselves tensors on spacetime and would 
therefore be integrated over spacetime using the measure 

(V 3 ?^ 0123 )- 

A comment on the dimensions of the variables. It 
would seem natural to measure the mass to and energy- 
momentum p of individual particles in the same units 
one measure gravitational mass and energy-momentum 
in, for example in the stress-energy tensor. But parti¬ 
cle momentum appears in the stress-energy momentum 
tensor only under the integration over momentum space. 
Dimensional analysis is therefore less restrictive than one 
might assume. 

It is consistent to assign particle mass and energy- 
momentum a dimension P that is independent of the 
dimension M of gravitational mass. In units in which 
c = G = 1, the unit of gravitational mass M is equal to 
the unit of length L, but P remains independent. For 
example, to. and p M have dimension P, and angular mo¬ 
mentum squared F = \x Apj 2 has dimensions L 2 P 2 , but 
the stress-energy tensor has dimension L , which can 
be thought of as ML~ 3 . The distribution function re¬ 
lates both kinds of masses and has dimensions L~ 2 P _4 . 
The particle current has dimension P ~ 1 L~ 2 , which can 
be thought of as MP~ l L ~ 3 : it has the dimension L ~ 3 of 
a number density only if gravitational and particle mass 
are measured in the same units. 

B. Spherical symmetry 


d dx»d dp i d i „ aJL 

da da dx ^ da dp 1 ^ dx ^ ^ dp 1 

(!) 

In this article a is an affine parameter related to proper 


Now we impose spherical symmetry. We describe the 
metric using polar-radial coordinates, 

ds 2 = —a 2 {t,r)dt 2 +a 2 (t 1 r)dr 2 +r 2 (dd 2 +sm 2 9d(j) 2 ). (5) 
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In this coordinate choice there is still the gauge freedom 
t —> t'(t), which changes the lapse a. We fix this freedom 
by setting a = 1 either at r = 0 or at r = oo. 

The Einstein equations give the following equations for 
the first derivatives of the metric coefficients: 


OL^ r 

a 

a t r 

a 

a,t 

a 


V + r r 8 ’ T - 

a 2 - 1 ra2 9, nT t 

- Y 8nTt ’ 


^8 

2 


( 6 ) 

( 7 ) 

( 8 ) 


The fourth Einstein equation, involving Tgg , is a combi¬ 
nation of derivatives of these three. The third of these 
equations, for a,*, becomes an identity when the other 
two equations are obeyed at each t. 

In spherical symmetry, the distribution function has 
to be of the form f{t,r,p r ,\p\), where \p\ 2 = gijp'p*. 
We can simplify the Vlasov equation using the fact that 
angular momentum is a constant of motion. Following 
H we define the variables 

w = ap r , (9) 

F = r 2 (\p\ 2 — w 2 ) = r 4 (p e “ + sin 2 9p^ 2 ). (10) 

We also introduce the shorthand 

W = \/m 2 + w 2 + F /r 2 = ap 4 . (11) 


In these variables the Vlasov equation in a spherically 
symmetric spacetime becomes 

df a w df f Fa a r a A df _ 

dt a W dr \r 3 aW a W a ) dw 

_ ( 12 ) 

where / is now a function f(t, r, w,F)- There is no 
derivative with respect to F because it is a constant of 
motion. 

The nonvanishing components of the particle number 
current are 


N* = 

. a- rv 

1-00 

/ dw /, 

(13) 


r z a J o j 

'-oo 

N r = 

= 4- f°° dF 
r 2 a J o J 

i- OO 

L dwf w ■ 

(14) 


The distribution function is always nonnegative and 
therefore we have aN l > a |lV r |, which means that N M 
is always timelike or lightlike. It is lightlike only for dis¬ 
tributions where all particles are massless and moving 
radially (F = 0) in the same direction, so that either 
w = W or w = —W. The non-vanishing components of 
the stress-energy tensor are 

nOO f- OO 2 

T/ = %l dF JA wf ^- 0 ’ (15) 

nOO p OO 

T/' = —- dF dw fW< 0, (16) 

^ JO J — oo 


r oo 1-00 

Tt = — 2 _ / dF dw f w, (17) 
^ ^ J 0 J — oo 

i-oo i-oo rp 

T ** = T ‘‘ = dF L dWf W^- (18) 

This tensor is conserved and satisfies the dominant and 
strong energy conditions. We assume that / behaves in 
such a way that all integrals (Hi) converge and are 
finite at every point. 

Assuming a > 1 the following inequalities follow from 
the Einstein equations: 


a r > 0, (19) 

(«a),r > |(« 2 ),t| , (20) 



Finally we have that 

i-OC 1-00 r 

= 0 . ( 22 , 


and therefore the Ricci scalar R is always nonnegative. 
For plotting results we shall also use the Hawking (or 
Misner-Sharp) mass function which is defined 

by 


a 2 (t, r) = 1 — 


2 M(t,r) 
r 


(23) 


C. Further symmetries and conserved quantities 

The Vlasov equation (J2j) expresses the fact that / is 
constant along trajectories of particles. Therefore the so¬ 
lution of that equation is, formally, an arbitrary function 
of the constants of motion of the problem. In the general 
case / would be a function of eight constants of motion, 
but when we impose spherical symmetry / can only de¬ 
pend on some of these. Spherical symmetry divides the 
four degrees of freedom of a point particle into two radial 
and two angular degrees of freedom. The angular part 
gives rise to four constants of motion: the three compo¬ 
nents L x , Ly, L z of the angular momentum vector, and 
the initial angle with respect to the axis given by this 
vector. If the averaged stress-energy of the particles is 
to be spherically symmetric, however, the distribution 
function / can only depend on the modulus (squared) F 
of the angular momentum. As we saw in the previous 
subsection, this can be used to reduce the Vlasov and 
Einstein equations to spherical symmetry. In the same 
way, additional symmetries can be used to further sim¬ 
plify the equations. 

The reduced, radial, system has two degrees of freedom 
[t(cr), r{<7 )] with conjugate momenta [pt(cr),p r (a)\ : where 
cr is an affine parameter along the particle trajectories. 
The Hamiltonian of the reduced system is 



( 24 ) 
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where F is now a given constant. This Hamiltonian sys¬ 
tem has four independent constants of motion c\, ..., 04 . 

One of these is the value of the Hamiltonian itself, 
2 

H = — The general solution of the equations of mo¬ 
tion can formally be written as 


t(a) = t(a,c 1 ,...,c 4 ,F ) 
r(cr) = r(a, Ci,c 4 , F) 
PM) = p t (a,c 1 ,...,c 4 ,F) 
p r (cr) = p r (a,Ci,...,C4,F) 


I => f = /(ci,...,c 4 ,F), 
(25) 


where we now consider the constants c.; as functions of 
t, r, pt, p r and er. Let c 4 = m. In the following we do 
not refer to to as a conserved quantity, but consider it a 
parameter of the equations. At least one of the four con¬ 
stants of motion, say C3, must depend on a, but because 
/ does not depend explicitly on a, f can not depend on 
this constant of motion. Therefore, the general solution 
/ in spherical symmetry will be an arbitrary function of 
F and two nontrivial constants of motion c\ and c 2 . 

Although the two nontrivial constants of motion c\ 
and C 2 always exist formally (and we can find them by 
integrating the equations of motion numerically), find¬ 
ing their analytical expressions in terms of the variables 
of the problem is only possible when the reduced radial 
system is integrable. Using the Liouville theorem [ flo| . 
we need to find only one constant of motion in terms of 
( t,r,pt,p r ), because then we will have two constants in 
involution for a system with two degrees of freedom. The 
easiest way of finding such a constant is imposing an addi¬ 
tional continuous symmetry. In the following we analyze 
two possible symmetries: staticity and self-similarity. We 
will then be able to give analytic expressions (at least 
quadratures) for the trajectories of the particles, sim¬ 
plifying the calculation of the energy-momentum tensor 
components. This is equivalent to giving a formal solu¬ 
tion of the Vlasov equation in terms of its characteristics 
in phase space. 

Even though we are interested in self-similarity, we first 
review the static case to show some important ideas in a 
well-known context. 


III. REVIEW OF STATIC SPACETIMES 
A. Field equations 

A static spacetime has an additional Killing vector 

S = 9t, (26) 

so that the metric functions a and a are just functions of 
r. The particles have an additional constant of motion 

E = —^Pfj, = —Pt = aW = ay/m 2 + w 2 + F / r 2 . (27) 

The four constants of motion (of the full 4-dimensional 
system) p 2 , F, L z and E are in involution and therefore, 
by the Liouville theorem, we can integrate the system. 


The radial equations of motion in terms of the constants 
of motion E and F are 


dr 

da 

dt 

da 


p r = ± 


E 2 


a(r) V a(r/ 


— TO 


E 

a 2 (r) 


F 
^2 > 


(28) 

(29) 


Defining v = p r /p* we can eliminate a and integrate the 
equations to obtain the quadrature 


where 


t-t 0 


r(t) dl j 


v{r') ’ 


(30) 




a(r) 

a(r) 



a 2 (r)F 

r 2 E 2 


We have the new constant of motion 


t 0 (t,r,E,F) = t- 



dr' 

v(r') ’ 


(31) 


(32) 


which is the time when the particle will be at ro, given 
that it is at r at time t. After E, this is the second of 
the nontrivial constants of motion Ci and C 2 referred to 
in the general discussion. 

Therefore, the most general solution to the Vlasov 
equation on a static background can be given as an arbi¬ 
trary function of three variables 


Vlasov in static spacetime => f(t, r, w,F) = g (to, E , F). 

(33) 

As to depends explicitly on t, this allows for a solution / 
that is explicitly time-dependent. 

If particles with a given E and F have bound orbits, 
g must be periodic in to with the period of that orbit. 

The period is determined as the integral (^0|) between 
the two turning points, and depends on E and F. We 
can avoid this complicated periodicity of g(to,E,F) by 
inverting equation (|3^) to obtain a constant of motion 
ro(t, r, E, F ) which gives the position ro of a particle with 
E, F at some canonical time to given that it is at r at 
time t: 


Vlasov in static spacetime => f(t, r, w,F) = g (ro, E, F ). 

(34) 

This means that we can freely specify f(r,w,F) at t = 
to, and still obtain an automatic solution of the Vlasov 
equation. 

So far we have considered Vlasov matter moving on a 
given spacetime. But to be consistent with staticity, we 
need the stress-energy momentum to be independent of 
t. Therefore, in principle, / should not depend on to or 
ro, which depend on t. The most general self-consistent 
static solution of Einstein-Vlasov in spherical symmetry 
would then have the form 


Einstein-Vlasov+staticity “ => ” f(t,r,w,F) = h(E,F). 

(35) 
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This result is referred to as Jeans’ theorem and it is 
known to hold in Newtonian physics (that is, in the 
Vlasov-Poisson system), but there are some counterex¬ 
amples in the Einstein-Vlasov system. That is why we 
have put the implication arrow in quotes. In order to un¬ 
derstand this subtlety we first give some explicit solutions 
where that result is valid. 


B. Some explicit solutions 


for the boundary and integral kernel. E m i n (r,F) gives 
the minimum energy of a particle at r with angular mo¬ 
mentum F, which corresponds to the state with w = 0. 
Note that it is always positive. It also functions as an 
effective potential for the radial particle motion, as we 
have 

v 2 = ^( E 2 ~ E L n )- (42) 


We choose a function h of two variables. The space- 
time can then be determined by solving the Einstein 
equations, which become a system of integral-differential 
equations. The integrals occur in forming the stress- 
energy tensor, which depends on the metric as well as 
on the arbitrary function h(E, F). 

In a static spherically symmetric solution we have two 
nontrivial Einstein equations. They are just the Einstein 
equations (jh)) and (|t|) in the general case: 


a 

a 


a 

a 


a 2 — 1 Air 2 ci 2 


2 r 


w 2 

dF / dwh(aW,F)—, (36) 

J — OO 


J 0 J — oo 

a 2 — 1 47r 2 a 2 

-1-x 

2 r r 


pOO poo 

dF dw h(aW, F) W, (37) 
J 0 J — oo 


where a and a are now functions of r only. In order to 
carry out the integrations, it is useful to change the inte¬ 
gration variable from w and F to E and F, the arguments 
of the free function h: 


a 

a 


a 

a 


2 r ,2 


a 2 — 1 47r 2 a 

+ 


2r 


p OO p OO 

/ dF 2EdE h(E,F) K 1/2 (38) 

Jo J Emintr.F) 


Emin(r,F) 

a 2 — 1 47r 2 a 2 


2 r 


ra- 


p OO pOO 

/ dF 2EdE h(E,F) K~ 1/2 , (39) 

Jo J Em\i\(r,F) 


' E mia (r,F) 

where we have used the shorthand expressions 


-Emin (r, F ) 
K(r, a, E, F) 



f O F \ 1/2 


a(r) 

1 ^ j > 

(40) 

w 2 

m 2 a 2 (r) 

Fa 2 (r) 

W 2 

E 2 

r 2 E 2 


E 2 - 


1 - 

E 2 

(41) 


When the function h(E, F) has been specified, the 
stress-energy integrals become given functions of r and 
a(r). The Einstein equations then become two coupled 
ordinary differential equations (ODEs) that determine 
a(r) and a{r). Conversely we can think of a(r) and a(r) 
as given functions. We then have an integral equation 
for h(E, F). By function counting, this integral equation 
does not determine h(E , F ) uniquely. There is an infi¬ 
nite number of functions h(E, F) that give rise to the 
same static spacetime, with the same stress-energy ten¬ 
sor. The change of integration variables to E and F 
shows that, for given h(E,F), the integrals in (pq) and 
( p9| ) are functions only of the two combinations r/a and 
ma. In particular, when m = 0, those integrals for fixed 
h(E, F ) are functions only of the combination r/a. 

The particular class of functions h(E , F) where 
h{E,F) = d(E 0 - E)9(F - F 0 )^(E)(F - F 0 ) 1 (called a 
“polytrope”), with upper cutoff E 0 and lower cutoff Fq, 
and l > —1/2, is analyzed in detail in [ [l~L| for massive par¬ 
ticles, where it is shown that solutions exist that describe 
matter distributions with finite total mass and compact 
radial support. 

If h(E, F) is restricted to a function of E alone, one 
finds by direct comparison of the integrals in the stress- 
energy tensor that T r r = Tg e = T ip ‘ p , so that the pressure 
is isotropic. This means that the stress-energy content in 
this case (spherically symmetric and static) has perfect 
fluid form. If p = —Tj and p = T r r are monotonic, one 
could also read off a, fairly meaningless, formal “equation 
of state”. 

It is interesting to solve a simple example in closed 
form, namely the ansatz h(E,F) = c9{E — Eq), with c 
and Eo positive constants. This is a “polytrope” with l = 
0 and no cutoff in angular momentum. For convenience 
of notation we define the dimensionless constant fh = 
m/E 0 , and the constant c = cEq. With the gauge choice 
a(r = 0) = 1 we have E > m, and therefore m must be 
in the range from 0 to 1. The parameter c can always 
be set to 1 by a choice of length units, and is therefore 
trivial. The integrals in equations ( ]38| ) and (^9|) can be 
evaluated explicitly, and give 
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P 

V 


—Tt = 

T i r _ 

r — 


7TC 
2 ^ 
7rc 


\/l — m 2 a 2 (2 — m 2 a 2 ) + m 4 a 4 log 


1 + i/l — m 2 a 2 


2a 4 


— \/l-mV(2 — 5m 2 a 2 ) — m 4 a 4 log 


1 + -\/l — m 2 a 2 


(43) 

(44) 


for ma < 1, and p = p = 0 otherwise. For massless 
particles, these expressions simplify to 

p = (45) 

or 

so that we have the formal equation of state p = p/3. 
a(r) is an increasing function, defining the surface of the 
“star” at ma = 1. Note that lighter particles give larger 
stars. 

An explicit numerical example for the resulting metric 
is shown in Fig. |I], both for massive (m = 1/2, continuous 
line) and massless (m = 0, dashed line) particles. As 
we expected, the massive case gives a star of finite mass 
and finite radius, while the massless case gives an infinite 
mass, infinite size distribution of matter. This can be 
easily understood in terms of the motion of the particles: 
in this example massive particles follow bounded orbits, 
so that we have a static situation with a finite number 
of particles in the system, but massless particles follow 
unbounded orbits, so that we need an infinite number of 
them to get a static distribution. 

C. Counterexamples to Jeans’ Theorem 

We shall now review some static solutions that do, af¬ 
ter all, depend on the third constant of motion tq. As¬ 
sume a static distribution of matter and its corresponding 
static gravitational field. Two different test particles can 
have the same energy E and angular momentum F at 
different values of r. It has been proved that in Newto¬ 
nian physics (the Vlasov-Poisson system) the distribution 
function must be the same at both points for those values 
of E and F |l^|. This result is referred to as Jeans’ the¬ 
orem. We can understand it physically in the following 
way. In Newtonian physics it is always possible to find 
a trajectory for one of the particles on which it reaches 
the position of the other particle without changing E or 
F. Because / is constant along particle trajectories it 
must be the same at both positions. This result is not 
true in general relativity: it is possible to show that two 
particles with the same E and F can be separated by a 
potential barrier. / can then be consistently different at 
those positions, and therefore it is not just a function of 
E and F, but of ro as well. 

The key observation is that the Newtonian minimum 
energy 

£min ( r,F) = U (r) + ^2 > ( 46 ) 


where the gravitational potential U (r) obeys the Poisson 
equation, always has a single minimum and no maxima, 
while the relativistic function may have several extrema. 

Rein has given a nice example in [jlH : Given a 
Schwarzschild black hole of mass M, it is impossible to 
set up a static, finite-mass distribution of Vlasov particles 
very close to it (for r < 3 M... 6 M, depending on their 
angular momentum) because it is not possible to have 
bound stable orbits in that region. Therefore we must 
set h = 0 in that region. However we can have a static 
distribution of particles with the same energy and angu¬ 
lar momentum further away from the black hole. The 
two regions are separated by a potential barrier. There¬ 
fore h may be different from zero for the same values of 
E and F. Again, h depends on r 0 , as well as on E and 
F. The left side of Fig. || gives an explicit numerical 
example of a self-gravitating shell of collisionless matter 
surrounding a black hole. 

A second example has been given by Schaeffer (li| . He 
shows that a shell of matter generates a well in the po¬ 
tential which can sustain another shell, if the inner shell 
is concentrated near its own Schwarzschild radius. We 
have an example in Fig. || on the right. The second shell 
can have a different value of the distribution function h 
for particles with the same energy and angular momen¬ 
tum. In the example in Fig. |] h in the second shell is half 
of h in the first shell. The same values of E and F are 
possible in a third region stretching to infinity (unbound 
particles), and h has been chosen to vanish there. 

We now state the Einstein-Vlasov form of “Jeans’ the¬ 
orem” . In general relativity particles with given values 
of E and F can exist in disjoint regions n = 1,2,... for 
the particles. Each region is bounded by turning points 
t±\e,F). We assume that r^ 4 < r < // +1 ' 1 . If the first 
allowed region contains the center then r^ = 0 and if 

the last allowed region is not bounded then 189,4 = oo. 
The general form of the distribution is 

h(r 0 , E,F)=J2 %o - rL n) ) - r 0 ) h^(E, F). 

n 

(47) 

In order to find self-gravitating solutions of this form, 
we shall usually need a convergent iterative procedure. 
For example the shells in the previous counterexample 
given by Schaeffer can be constructed iteratively using 
the fact that in spherical symmetry the interior shells are 
not affected by the exterior shells (Birkhoff theorem). 

The dependence of h on ro in ([47]) is admittedly rather 
restricted, but we have discussed it in some detail, be- 
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M' a 




FIG. 1: Mass distribution M'(r) (left) and lapse function a(r) (right) corresponding to the solution given in Eqs. ([1^) and 
(0 for m = 1/2 (continuous line) and m = 0 (dashed line). 


E m \n/lTL E m \ n /7Tl 




FIG. 2: Graphs of the effective potential £'min(r, Fo) = a(r)\Jm 2 + Fo/r 2 , where Fo is a lower cutoff in angular momentum. 
The horizontal dashed line gives the upper cutoff Eo in energy. We can have particles between both lines (shaded regions). In 
both examples, m = 1, and h is independent of E and F for F > Fo and E < Eo, but takes different values in different potential 
wells. In both examples we have set a = 1 at r = oo, and so F? m in/m —> 1 from below as r —> oo. In the first graph we have 
a black hole with M = 1 and we plot just the exterior region. Fo = 17. We could have particles in the region 2M < r < 3 M, 
but not a static distribution of them. In the second graph we have two shells of matter, and Fo = 2. Both of them contain 
particles, but h is half in the outer well of what it is in the inner well. 


cause the dependence of the equivalent of h in self-similar 
solutions on all three constants of motion will be non¬ 
trivial, and in fact no solutions can be found that depend 
only on F and the equivalent of E: solutions must also 
depend on the equivalent of rg. 


D. Massless particles 

Massless particles are not very interesting when one 
looks for static solutions because it is difficult to form 
asymptotically flat solutions, but the system with mass¬ 
less particles has an additional symmetry: the trajec¬ 
tory of a massless particle, in a given spacetime, for a 
given initial position, is completely determined by the 
direction of its initial four-momentum, independently of 


the modulus. In order to make a given contribution to 
the stress-energy tensor, we can therefore use N particles 
of four-momentum p M or one particle of four-momentum 
Np^. We shall see that this allows us to state explicitly 
which part of the Vlasov function h(E, F) is determined 
by the spacetime metric and which is arbitrary. This is 
interesting as a toy model for the CSS solutions (where 
we are interested in massless particles anyway). 

In a spherically symmetric static solution in particular, 
the trajectory of a particle is determined not by E and 
F separately, but only by the combination E 2 /F, which 
is invariant under rescalings of the four-momentum. We 
rewrite the free function h(E, F ) as 

h(E,F) = k(y,F), 


(48) 



















y = 


E aW (w 2 1 \ 2 

7f = 7f = a \f + 


(49) 


With m = 0, the limit E m ; n (r, F) in the integration over 
E becomes independent of F , and the two integrations 
can be interchanged. We find 

a' a 2 — 1 47r 2 a 2 

— = —o-1- 2 x 

a Zr ra z 

f Hy) 2 v d y > (50) 

J a /r \ r z y z ) 

a! a 2 - 1 47r 2 a 2 

— = o-'-2 x 

a Zr ra z 

[ k(y)2ydy, (51) 

J a /r\ f y J 

where k is the integral of the Vlasov function k over all 
particles with the same trajectory: 

/*oo 

k(y)= / k{y, F)F dF. (52) 

Jo 


The integration limit y = a(r)/r is the turning point 
(perihelion) of all particles with a given y. Note that 
[/] = Ill] = [h] = [&] = L~ 2 P~ 4 (for both massive and 
massless particles), and [k] = L 2 . 

For massless particles the Ricci scalar vanishes, and as 
a consequence a(r) and a(r) are not independent but are 
related by a second order ODE. This ODE has a unique 
solution for either a(r) given a(r ) and the boundary con¬ 
ditions a(r) = 1, a'(r) = 0, or the other way around. a(r) 
and a(r) are therefore related one-to one. 

Furthermore, Eqs. (^0|M) can be turned into linear 
Fredholm equations of the second kind by a change of 
variable. This means that either of these equations can 
be solved for k{y) for given a(r) and a(r). The relation¬ 
ship between a(r) and a(r ) on the one hand and k(y) on 
the other is therefore one-to-one. 


IV. CSS SPACETIMES 
A. Spacetime geometry 

Because the particle mass m introduces a scale into 
the field equations, it is not clear that massive parti¬ 
cles are compatible with exact self-similarity. In this pa¬ 
per we therefore construct CSS solutions of the massless 
Einstein-Vlasov system. We defer the question if mas¬ 
sive particles allow CSS solutions or asymptotically CSS 
solutions to another paper. 

We start by defining CSS and stating some general 
properties of CSS spacetimes. Then we analyze the possi¬ 
ble trajectories of massless particles in those spacetimes. 
We find two nontrivial conserved quantities for massless 
particles, J and tq , which are the equivalents of E and to 
in the static case. In contrast to the static case, we will 


find that the Vlasov distribution function / can not just 
depend on J and F but must also depend on To in order 
to construct a regular self-similar stress-energy tensor. 

A CSS spacetime is one that possesses a homothetic 
vector field 


+ V,,^ — —2 g^. (53) 

In spherical symmetry, in the coordinates r and t, this 
vector field is 


£ = ~td t - rd r . (54) 

£ is the infinitesimal generator of scale transformations 
r —> sr, t —> st. It is customary to define the following 
coordinates, which are adapted to self-similarity in the 
region (r > 0,t < 0): 

XS -f TS - bg(V) , (55) 

where Lq is a constant put in to get dimensions right. It is 
helpful to assign t, r, and Lq dimension L, while x and r 
are dimensionless. Fig. || gives a graphical description of 
the new coordinates. In these coordinates the homothetic 
vector is 


f = d T , 


(56) 


and the metric functions a and a are functions of x only. 
(We could write the metric explicitly in terms of coordi¬ 
nates x and r, but it is often more helpful to refer back 
to coordinates t and r as they are more familiar.) 

Geometrically, the self-similar spacetime singles out 
the point r = t = 0, which in fact is a curvature singu¬ 
larity. It will be useful to introduce two related auxiliary 
functions, 


G(x) 


a(x) x 
a(x) 


H(x) = USE! - a 2 {x). (57) 

x z 


The past light cone of the singularity (from now on also 
referred to as “the” light cone) is the set of points (r, xi c ) 
with a^ic the solution to the equation 


G{x\c) = 1 . 


(58) 


(Of course G = 1 is equivalent to H = 0.) In the fol¬ 
lowing, we assume that there is only a single light cone 
G = 1 in the region t < 0. 

Once we know the behavior of the metric under scale 
transformations we can derive the behavior of the Ein¬ 
stein tensor, and therefore of the stress-energy tensor. 
The components of the stress-energy tensor in coordi¬ 
nates (t,r) obey 

T/(f,r) =t- 2 T/(x). (59) 

Using the rescaled components TjA, the Einstein equa¬ 
tions become ODEs in x : 


a 

a 


a 2 - 1 


2 x 


47 rxa 2 Tt (x ), 


(60) 
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FIG. 3: Map between coordinates (t, r) and (t,x). The past light cone of the point t = r = 0 has been arbitrarily set to 
x = xic = 2. 


-1-= 47rxa 2 T r r ’(a;), (61) 

a 2x 

— = ATTa 2 T t r {x). (62) 

a 

Because in CSS a is a function of x only, the equations 
for a : t and cq r have both become equations for a', and 
can be combined to eliminate a' and obtain an algebraic 
equation for a, 

a~ 2 (x ) — 1 = 87t \x 2 Tt{x) + xT t r (x)\ (63) 

In a CSS spacetime, the Hawking mass M ( t , r) is 
of the form M(t,r) = e~ T M(x). At the regular cen¬ 
ter, M{ 0) = 0, and from the Einstein equations it fol¬ 
lows that M'(x) > 0. Suppose now that a spheri¬ 
cally symmetric CSS spacetime contained a vacuum re¬ 
gion, so that no matter was present in some interval of 
x. The only spherically symmetric vacuum solution is 
Schwarzschild. But a Schwarzschild solution has a finite 
Hawking mass M(t,r) = M, which is incompatible with 
M(t,r ) = e~ T M(x), unless M(t,r ) = 0. This means 
that a spherically symmetric regular CSS solution can 
contain a central vacuum region, but no vacuum region 
surrounding a matter region. In particular, it can not 
be asymptotically flat. (The exception would be a space- 
time with a negative mass naked singularity at the cen¬ 
ter, as M(0) < 0 would make M(oo) = 0 compatible with 
M’(x) > 0.) 

We note in passing that for massless particles the 
energy-momentum tensor is traceless, and therefore both 
the Einstein and Ricci tensors are traceless. This gives 
an additional relation between the metric functions a and 
a: 



Thus, we only have to solve for one of the metric func¬ 
tions, or equivalently, we only have to know one of the 
components of the energy-momentum tensor. 


B. Massless and massive particles 


We now discuss test particle trajectories in a fixed CSS 
spacetime in some detail, as their properties are unfamil¬ 
iar. Starting from the homothetic vector (^3|) is natural 
to construct a “homothetic energy”, 


J = = -p T = tp t + rp r . (65) 


In a self-similar spacetime it obeys the equation of motion 


dJ_ 

da 


= —m 


2 


( 66 ) 


where a is the affine parameter along particle trajecto¬ 
ries. For massless particles J is a constant of motion. 
For massive particles, J + m 2 a is constant, but to use 
this constant of motion one would have to solve for a as 
a function of (t, r, w, F). 

In a CSS background, it is natural to rewrite the Vlasov 
equation in terms of f(t, r, J, F), obtaining 


t df df 2 df 

v m +r, S-r- m gj = 0 ' 


(67) 


where p r and p t are now functions of (f, r, J, F). 

In order to integrate trajectories of the reduced system 
in r and t, we define 


dr _ p r 
dt p* 


( 68 ) 


p r and pt can be obtained from F and J by solving the 
system 

J = —a 2 tp t + a 2 rp r , (69) 

F = r 2 (a 2 p t2 — a 2 p r2 — m 2 ^j . (70) 

We know that F > 0, and by convention p 4 > 0. The 
radial momentum p r can be positive, negative or zero, 
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and therefore J may also have either sign. For a(x)x > 
a(x), that is outside the light cone, a unique solution 
exists for all J. For a(x)x < a(x), that is inside the light 
cone, one solution exists for J = J c , and two solutions 
exist for J > J c , where 

J c (x, F,mr) = y/H(x)(F + m 2 r 2 ). (71) 

The condition J > J c means that we are outside or at 
the minimum radius that a particle with given x and 
F can have at a given t. On the light cone, a(x)x = 
a(x), we have a unique solution for J > J c = 0. On 
the light cone J = 0 is only possible with F = m = 0 
and applies to massless particles of arbitrary momentum 
moving radially into the singularity along its past light 
cone. 


C. Massless particle trajectories 


We now discuss the trajectories of massless particles in 
more detail. J is then a constant of the motion. From 
the definitions of v, J and F it is clear that for m = 0, v 
depends only on x and the combination 

rm W <72) 

With the convention p t > 0 we find that 

y = £ yx+£L <73) 

\/a 2 — a 2 v 2 

To solve this equation for v(x,Y), we need to square it 
first. This gives rise to a sign ambiguity that we discuss 
below. The result is 


v = 



—aa ± Y\J—a 2 + (a 2 + Y 2 )x 2 
a 2 + Y 2 


v±(x,Y). 


(74) 

The allowed range of v for particles with F > 0 is—a/a < 
v < a/a. The rates of change of x and r are related by 


— = V = v + x (75) 

dr 

In coordinates x and r, the particle equation of motion 
is autonomous: 

= v±(x,Y) +x = V±{x,Y). (76) 

dr 

Trajectories of massless particles are therefore deter¬ 
mined by y, up to a translation in r. A translation in 
r corresponds to a simultaneous rescaling of t and r. A 
second symmetry of massless particle trajectories arises 
because Y is invariant under a simultaneous rescaling 
of all components of This leaves the trajectory in 
spacetime completely invariant: photons of different en¬ 
ergy can have the same trajectory. The overall momen¬ 
tum scale, at constant Y and constant r and t, is set 


by y/F. Here we assume that the contribution of parti¬ 
cles with F = 0 to the stress-energy is zero, and we now 
concentrate on particles with F > 0. (If they made a fi¬ 
nite contribution, this would constitute null dust. In our 
numerical examples we avoid this simply by imposing a 
cutoff.) 

The projection into the (r, t) plane of a the trajectory 
of a massless particle with F > 0 is timelike. Such par¬ 
ticles come in from infinity, reach a minimum radius r, 
and go out to infinity again. At large r, the tangential 
velocity can be neglected compared to the radial velocity, 
and the radial projection of the trajectory is asymptoti¬ 
cally null. Fig. H shows several trajectories in the (r, t) 
plane that share the same value of Y but are related by 
a rescaling of r and t. In the (a;, r) plane these would be 
related by a translation in r. Fig. |5j shows typical par¬ 
ticle trajectories in the (r, t) plane with different values 
of Y. Fig. shows the same trajectories in (x, r) coor¬ 
dinates, and Fig. [ 7 ] shows the equivalent trajectories in 
(x, V) space. One should keep in mind that the turning 
point of smallest x (if one exists) does not coincide with 
the turning point of smallest r. 

Fig. H illustrates that there are two types of particles: 
type I particles are initially inside the light cone but leave 
it, while type II particles are always outside the light 
cone. Both types of particles start at x = x\ c in the 
asymptotic past as r —> — 00 . Type I particles peel off 
from x = X\ c towards smaller x, reach a minimum value 
x m in (Y ) of x and turn round. A finite r-time later they 
cross the light cone x = x\ c , and continue to large x. 
Type II particles peel off from x = x\ c towards larger 
x. Although they reach a minimum value of r, their 
value of x always increase. Both types of particles reach 
x = 00 at r = 00 . This is just the surface t = 0, where 
our coordinate system breaks down. Here we are only 
interested in the region t < 0. 

In order to see when dx/dr vanishes, and which of the 
two values V± is the correct one in a given situation, it 
is helpful to factorize V± as 


V±(x,Y) = x(a 2 +Y 2 )~ 1 (Y 2 —H ) 1 ^ 2 (Y 2 — if) 1 / 2 ± G~ X Y 

(77) 

where G(x) and H(x) were defined in (£7|). Clearly both 
V+ and V- vanish at x = x m [ n (Y) which is defined by 


H(x mia ) = Y 2 


(78) 


One can show that at these points d 2 x/dr 2 does not van¬ 
ish and has the same sign as — dH/dx. Therefore these 
points are true turning points. The square bracket in 
( f77l ) can vanish only at the light cone x = xi c , where 
H = 0 and G = 1. In particular V+ vanishes at x = x\ c 
for all Y <0 (but not for Y > 0), and V_ vanishes for all 
Y > 0 (but not for Y < 0). This corresponds to V —> 0 
as r —> — 00 (compare Fig. |(ij), and one can show that 
d 2 x/dr 2 = 0 there. 

To see how the sign choice V± relates to our classifica¬ 
tion of particles, we note first that Y and — Y give rise to 
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t 




FIG. 6: The same trajectories in coordinates (x,t). From 
left to right we have Y = 10, 5, 3, 2,1, 0, —1. The thicker line 
is Y = 0. The origin in the r axis has been arbitrarily chosen 
(corresponding to a choice of units Lq). 


FIG. 4: Trajector ies of particles in the self-similar spacetime 
described in Sect. IV E| below, all with the same value Y = 5 
but different overall scales. The past light cone of the sin¬ 
gularity at xic = 1.3642 is denoted with a short-dashed line, 
and for orientation x = 2 has been marked with a long-dashed 
line. 


t 



FIG. 5; _ (t, r)-trajectories of particles in the self-similar space- 

time IV E with Y = 10, 5, 3, 2,1, 0, —1 (from left to right). 
The trajectory with Y = 0 has been indicated by a thicker 
line. The overall scale of each trajectory has fixed so that all 
trajectories meet in one point. 



0 0.5 1 1.5 2 


x 


FIG. 7: The equivalent trajectories (from left to right Y = 
10, 5, 3, 2,1, 0, —1) in variables x and V. 


cone. 

From Fig. ^ we see that type II particles have v < 0, 
and therefore p r < 0. (To be precise, this holds for the 
region t < 0 covered by our coordinates. Presumably 
these particles have an r-turning point for some t > 0. 
Y = 0 particles have an r-turning point at t = 0.) There¬ 
fore type II particles have J < 0, and so F <0. This is 
summarized in Table ||. 

It seems plausible that CSS spacetimes exist where 
H{x) has a local minimum, so that for some range of Y 2 
the turning point equation H = Y 2 would have two so- 


inequivalent trajectories. With t < 0 and p * > 0, it fol¬ 
lows from ( ]69| ) that J > 0 for any particle inside the light 
cone. Therefore type I particles have Y > 0. Inside the 
light cone, the velocity of type I particles is V = VL < 0 
while they are ingoing, and V = V+ > 0 while they are 
outgoing. By continuity, V = V + also outside the light 


TABLE I: Type I and type II particles 



Type I 

Type II 


Y > 0 

Y < 0 

inside light cone 

V = V± 

never there 

outside light cone 

v = v+ 

II 
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lutions £ m in and a; max for each value of Y. Particles with 
Y in that range would oscillate between the two turn¬ 
ing points. Such “type III” particles would be trapped 
within a finite range of x, and would therefore hit the 
singularity. However, from the definition (|6^) of J it fol¬ 
lows that p r and p 4 at constant x and J scale as t _1 . As 
a type III particle returns to the same value of x again 
and again, its energy and radial momentum as measured 
by a constant r observer increase as t -1 . A similar argu¬ 
ment holds for the tangential momenta. A finite number 
of such particles would therefore constitute an energy- 
momentum scaling as f -1 confined to a region of a physi¬ 
cal size scaling as t , so that the energy and other frame 
components of the stress-energy would scale as t~ A . This 
is incompatible with self-similarity, where these stress- 
energy components scales as t ~ 2 . Type III particles are 
allowed as test particles, but can not be used as the mat¬ 
ter content of a CSS spacetime. From now on we restrict 
consideration to type I and II particles. 


D. Solutions with massless particles 


For massless particles, J is a constant of motion, and 
therefore d/dJ disappears from the Vlasov equation (67). 
We then obtain an automatic solution of the Vlasov equa¬ 
tion by the ansatz f(t,r,J,F) = h(J,F). However, we 
shall see that this ansatz leads to a divergent stress- 
energy tensor. Therefore we have to resort to the most 
general situation on a CSS background, by allowing for 
the dependence of / on the spacetime coordinates r and 
t through a dependence on a third non-trivial constant of 
the motion. (With ci = J, this is the constant C 2 of our 
general discussion of conserved quantities.) We rewrite 
the Vlasov equation for massless particles in coordinates 
x and r as 


O n 

— +V(x,Y) — ) f(r,x,J,F) = 0, (79) 


We can solve this by the method of characteristics as 


/( T > x,J,F) = g(r 0 ,J,F), 


where 


tq{t,x,Y) = t - 


dx' 


0 (X) V{x',YY 


(80) 


(81) 


Clearly r 0 is the time r when the particle with a given 
Y and F was at x = xg- Here we have specified xo as a 
given function of Y (which could be a constant). A priori 
xq could also depend on F. but that would break self¬ 
similarity of the resulting stress-energy. For xo (Y) < x\ c 
(only possible for type I particles) we must specify if the 
spacetime point (xo,Tg) is on the ingoing or outgoing 
part of the trajectory. A natural choice then is Xo(Y) = 
Xmin(Y). but this does not work for type II particles, 
which do not have a turning point. For either type of 
particle we could use a fixed value Xg = const > x\ c . 


As we are mainly working in coordinates t and r in this 
paper, we shall use r 0 in the form of 

r x dx' 

t 0 = -Lge T ° = t exp / = tQ(x,Y). 

Jx 0 (Y) V(x',Y) 

(82) 

(Note again that Q(x,Y) depends on the choice of 
x o(Y)-) 

We have just shown that the general spherically sym¬ 
metric solution of the Vlasov equation in a spherically 
symmetric CSS spacetime is 


Vlasov in CSS spacetime => f(t, r, w, F ) = g (to, J, F ). 

(83) 

One could think, by analogy with the static case, that 
the requirement of a self-similar energy-momentum ten¬ 
sor eliminates the dependence on to- However, that is 
not the case. In the static case, the constants of motion 
E and F depend only on the variables r, p* and p r . Only 
the constant of motion to depends also on t, and that 
is why g can not depend on it (except in the piecewise 
constant way of Eq. (|47|)). In the self-similar case, not 
just to but all three constants to, J and F change under 
a scale transformation. Imposing consistency with a CSS 
spacetime is therefore less straightforward. 

In order to derive the appropriate form of g(to, J, F), 
it is helpful to change the integration over F, which 
has dimension P 2 L 2 , temporarily to an integration over 
F = F/r 2 , which has dimension P 2 . This completely 
separates the integration over momentum space from the 
remaining dependence on r and t. The rescaled stress- 
energy tensor components are 


Tr 


r 


ft 


t 


Tt r 


T 9 ° 


7T j dF J dwt 2 g(t 0 ,J,F)^y 


—7 r / dF / dwt 2 g(to,J,F)W , 

J 0 J — oo 

^ p oo poo 

—7r— / dF / dwt 2 g(t 0 , J, F)w, 
a Jo J-oo 

p OO p OO pi 

2 L dP _J«t 2 9(t a ,J,F) w - 


(84) 

(85) 

( 86 ) 
(87) 


The functions to, J , F, and W are homogeneous func¬ 
tions of r and t of degree 1, 1, 2 and 0, respectively. If we 
now demand that T^ v (st,sr) = f tl v (t,r) for any scaling 
constant s, so that the stress-energy tensor is compatible 
with CSS, we find that g must be the following homoge¬ 
neous function of its arguments: 

Vlasov + CSS => g(sto, sJ, s 2 F) = s~ 2 g(tg, J, F) 

( 88 ) 

A homogeneous function of three variables can be re¬ 
expressed as an arbitrary function of two variables. We 
define 

9 (to, J, F) = -U’Qb z ), where Z = ^ (89) 

r £ 0 

and where Y was defined above in (f72|). The situation is 
similar to static solutions of Einstein-Vlasov in that we 
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can freely specify a function k(Y, Z) of two variables that 
automatically solves the Vlasov equation, and obtain the 
spacetime by solving the Einstein equations, which be¬ 
come integral-differential equations. It differs in that 
g(to,J,F) is not simply of the form h(J,F). Because 
we need to in order to obtain the general solution, we 
need to determine Q(x,Y), which is equivalent to inte¬ 
grating the equations of motion. Note that f(t,r,w,F) 
is determined by h(Y,Z) only when the function Xq (Y) 
that appears in Q has been specified. 

It is interesting to note that the homogeneity condi¬ 
tion ( |S8|) arises from dimensional analysis in the length 
dimension L alone. Y has neither P nor L dimension, 
but Z has dimension P 2 . This is not an obstruction to 
self-similarity of the spacetime. If we had measured par¬ 
ticle energy-momentum in units of length, dimensional 
analysis would have confused us in this point. Note that 
[/] = [5] = L~ 2 P~ 4 , [k] = P~ 2 , and [k] = L 2 . 

The Einstein equations are now once again a set of 
integral-differential equations that we can solve by iter¬ 
ation, once k(Y,Z) and Xo(Y) have been specified. The 
next step in our program would be to change the inte¬ 
gration variables from F and w to Y and Z. By analogy 
with the massless static case, we expect the integration 
over Z to decouple from that over Y. This is correct, 
but going to the variables Y and Z directly gives rise 
to rather complicated expressions. When calculating the 
stress-energy components, we find it helpful to integrate 
over the variables Z and u = w/\fF. This already de¬ 
couples the integrals: 

t r r°° _ „ 2 

Tr r = - / Q 2 {x, Y) k(Y)—j===du, (90) 

X 4 J- oo Vu 2 + 1 

roo _ 

Tt* = —7 Q 2 (x,Y)k(Y)Vu 2 + ld Ul (91) 
X 4 J-oc 

/ OO 

Q 2 (x, Y) k(Y)u du, (92) 

-00 

, 0 ° _ 1 
L q1(x ' y)Hy) 7^tV‘ u - (93) 

where 

/»oo 

k(Y)= / k(Y, Z) dZ. (94) 

Jo 

The auxiliary variable u is related to Y by 

Y = Y(x, u) = y/u 2 + 1 + a(x)u. (95) 

x 

We can now see clearly why g must depend on to- If it 
did not, k(Y. Z) could not depend on Z, and therefore 
k(Y) would diverge. The dependence on to is necessary 
to create a finite CSS stress-energy tensor. 


E. Some explicit solutions 


Q±{x,Y) can then be calculated in closed form. From 
the definition of Q it follows that it obeys the differential 
equation 


dlnQ±(x,Y) 1 

dx V±{x, Y) 


(96) 


Integration of this is straightforward when a(x) = a(x) = 
1. We now consider the boundary conditions. We make 
the assumption that the reference point ( Xq , To) is on the 
outgoing branch of the trajectories. This gives rise to the 
boundary condition Q + (xo(Y), Y) = 1 for outgoing par¬ 
ticles. For ingoing particles we have to piece Q together 
from an ingoing and an outgoing piece, integrating first 
from x to x m i n using V- and then from x m i n to Xq using 
V + , respectively. Here the turning point :r m i n of type I 
particles is given by 

^min(V) = (1 + V 2 ) -1 / 2 (97) 


for Y > 0. (Type II particles, with Y < 0, have no 
turning point.) Our final result is 


Q±(x,Y) 


Y±y/a?{l + Y 2 ) - 1 
Y+y/x%{Y)[l+Y 2 )-l' 


(98) 


which holds for type I and type II particles as long as x 
and xo{Y) are both in the physical range: x > £ m i n (Y) 
for Y > 0, and x > 1 for Y < 0. Note that Q_ —> 0 as 
x —> 1_ for type I particles and Q + — * 0 as x —> 1+ for 
type II particles: this is the limit in which the trajectory 
peels off the light cone at r —> — 00 . Note also that Q± is 
finite for type I particles as Xq —^► 1 because the reference 
point is the point where the trajectory crosses the light 
cone while going out. By contrast, Q + diverges for type 
II particles as xo —^> 1+ because this pushes the reference 
point to r —> —ex: where the trajectory peels off the light 
cone. 

On a flat spacetime, neglecting gravity, we are just 
dealing with the Vlasov equation, which is linear. It is 
therefore sufficient to consider the following ansatz for 
k(Y): 


k(Y) = k+ 6 {Y - Y 0 ) + k _ 6 {Y + Y 0 ), (99) 


with Yo > 0. We can then integrate the resulting stress- 
energy tensor over lo, with arbitrary functions k±{Yo), 
in order to obtain the general result. We shall see later 
why it is useful to consider the contributions Y = Yq and 
Y = —Y 0 together. 

We first calculate the stress-energy inside the light 
cone. Because type II particles are always outside the 
light cone, a solution with both type I and type II par¬ 
ticles is independent of its type II particle content for 
0 < x < 1. Inside the light cone we find that the stress- 
energy tensor and particle number current components 
are 

2t tK+ 1 

^min X\J X 2 ^min 


We begin with a CSS solution of the Vlasov system in 
fiat spacetime, which is trivially self-similar with x\ c = 1. 


Tt\x) = 


( 100 ) 
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Tr r (x) 


Tt r (x) 


T e e (x) 

N*(x) 

N r (x) 


2ttK + xJ n + (1 ~ x 2 min )(x 2 - J 2 lin ) 

^fnin X 3 ^/X 2 - X 2 min 

2ttK + 1 

Zmin SV* 2 -*mm’ 

nK + 1 + x 2 - 2x 2 min 

Xmin X 3 ^/x T ~xl~ , 

2ttK+ yjl - x 2 min 

X 2 ■ *t>2 /™2 _ ~2 ’ 

mm X x min 

2^+ V^-^min ^-^min 

^niin V* 2 — X min ‘ T 


( 101 ) 

( 102 ) 

(103) 

(104) 

(105) 


for x > £ m i n , and vanish for x < x m \ n . Here cc m j n is short¬ 
hand for ar m inQ'o)j and we have introduced the short¬ 
hands 


K± = 


k± 


(y 0 ±^x 2 (1 + Y 2 )-1 


2 > 


(106) 


These expressions already take into account the contri¬ 
bution of both ingoing and outgoing particles. We have 
used Q given by Eq. (^) so that the reference point 
(xo,to) is always assumed to be on the outgoing branch. 
All these expressions vanish for x < a; m i n , where no parti¬ 
cles reach, and have an integrable divergence at x = a; m i n , 
where particles pile up as they turn around. 


We now calculate the stress-energy outside the light 
cone. For type II particles, we must use xo(Y) > X\ c = 1. 
In the following, we use the same value Xq > 1 for both 
xo(ho) and Xo(—Yq). It is therefore automatically on the 
outgoing branch of type I particles. Taking into account 
the contributions of both type I and type II particles, we 
find for x > 1 that 


Tt\x) 


T r r {x) 


Tt r (x) 


T e e (x) 

N\x) 

N r {x) 


K+ + AT 

xsjx 2 - x 2 min ’ 

(K -A- K 1 ^ (1 WninX 3 -' a 'min) ■) ( 'j ^min 1 

+ x 3 Jx 2 -xY + 


X 


min 


(K+ + K-) J. - - (K + - K_) 

x 1 \Jx 1 -xY n 


y/1 ~ x 2 


min 


min 

(K+ + K-) 3 + /2 2 mm + 2 (K + - KJ) 

X W X2 -<in 

(k + +kj )4 + (k + -kj) = 

* X 2 V x2 ~~ 

1 _ p T 2 /i _ t 2 „2 _ o™2 

(K+ + KJ) - pJs. + (K+ - K ) v 

X \/x 2 - x'j X 


2 2 
^min 


(107) 

(108) 

(109) 

( 110 ) 
(HI) 
( 112 ) 


The constant AT was defined above. Here, as above, 
we use x m in as shorthand for the turning point a; m in (Yo) 
of the type I particles. (This number also appears in the 
contribution of the type II particles, but is then just used 
as a shorthand.) 

Based on these results, we can make two important 
general points. The first is that for general values of K± 
the stress-energy tensor is not continuous at the light 
cone. However, if and only if A+ = A_ the stress-energy 
is continuous at the light cone, and in fact analytic. On 
the other hand, the particle-current is always continuous, 
but never analytic. In the general flat-space case, where 
k± are promoted to functions of Y >0, the stress-energy 
is analytic if K + (Y) = K_(Y) for all Y. 

The second important point is that while we can freely 


specify both k(Y,Z) and xq(Y), the location of the ref¬ 
erence point xo(Y) is not really physical. We have seen 
that the stress-energy in our example does not depend 
on the two numbers Xq and k + (or kJ) separately but 
only on the combination K + (or KJ). Generally, the 
Vlasov density f(t,r,w,F) is completely determined by 
one free function of Y and Z. In flat spacetime, and as¬ 
suming that xq(Y) = £o(|V|), we have shown that this 
free function is 

K{Y,Z) - - - - r. (113) 

(|F| + Vx 2 (Y)(l + Y 2 )-l) 

In the self-gravitating case we cannot give this function 
in closed form. Therefore, in the self-gravitating case we 
must just “fix the gauge” by fixing xq{Y) arbitrarily. 
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x 


FIG. 8: Discontinuity in — Tt l for the 5-function case with 
only type I particles (continuous line). The d ots denote the 
same quantity for the Gaussian function (114) with the same 
k+/y/2 7r = 1CP 4 , Vo = 6 and a = 1. The agreement is very 
good. 



x 


FIG. 9: Stress-energy tensor corresponding to Gaussian 

(114) on flat spacetime. Note the jump at x = x\ c = 1 . 


We shall now smear the 5-function in Y, and then 
couple this distribution to gravity. In the following we 
restrict our ansatz to type I particles. Inside the light 
cone this is no restriction anyway. It allows us to use 
Xo{Y) = X' m i n (y), which simplifies the numerical calcu¬ 
lations. We also want to start with a matter distribu¬ 
tion whose stress-energy tensor in flat spacetime is dis¬ 
continuous, because we shall see that the gravitational 
back-reaction makes it continuous. As an example, the 
jump in —Tt is shown in Fig. |8| for the case Yo = 6, 
k+ = 10 _4 v / 27t and fc_ = 0, using xq(Y) = a" m i n (E). 

Still in flat spacetime, we smear the 5-function into a 
Gaussian: 


_ *-+ r -(Y-Y n ) 2 /2a 2 


k(Y) = — e 


(114) 


with a minimum cutoff just above Y = 0 and a maximum 
cutoff far from Yq. It is not possible to calculate the 
results analytically, so that we calculate them with a C 
code. The divergence of the stress-energy at x rn j n (Y) is 
then smoothed out by the integration over Y. and the 
stress-energy components are very smooth apart from at 
the light cone. The results for the stress-energy tensor 
are shown in Fig. |d| for the case a = 1 with Yo = 6 and 
k + = 10 -4 -\/27r as in the 5-function example. 

Once we couple the matter to gravity we must solve 
the integral-differential equations by iteration. That is, 
we start with some initial metric, say flat spacetime, and 
calculate the stress-energy tensor of a set of particles for 
a certain function k(Y). Then we integrate the met¬ 
ric corresponding to this stress-energy distribution using 
Eqs. @|I|) and calculate a new stress-energy tensor, 


and so on, until the process converges. In Fig. 10 we 
show the metric that results from coupling the Gaus¬ 
sian k of the previous example, with parameters Yo = 6, 
(7 = 1 and fco = 10~ 4 -\/27r, to gravity. The convergence of 
the numerical method is demonstrated in Fig. Oil. Fig. 
IL3 compares the radial pressure profile of the Gaussian 


ansatz in flat spacetime with the Gaussian coupled to 
gravity. Even though we started with a discontinuous 
stress-energy tensor on flat spacetime, we now find that 
the self-consistent stress-energy is now C° (but not C 1 ). 
This is illustrated in Fig. [l3|. As a consequence, the met¬ 
ric is C 1 . Another check of the procedure is given by the 
Eq.(||). 

With very low values of k + we can form solutions with 
a metric which is very close to flat spacetime. On the 
other hand, with larger values of h+ we get solutions 
which are very close to horizon formation. For example 
with fc + = 3T0~ 4 -\/27r, the maximum value of a is greater 
than 1.45 (giving 2 M/r ~ 0.53), which is even bigger 
than the maximum value of the same function in the 
Choptuik or Evans-Coleman spacetimes. 

We can understand analytically why the coupling to 
gravity makes the matter stress-energy more regular. For 
particles that just peel off the light cone, that is for x ~ 
X] c . and V = V- for type I particles and V = V+ for type 
II particles, we can expand V as 


V±(x, t\Y\) = c(x - xi c ) + 0(x - xi c ) 2 , 


c = x\ c G'{x\ c ) 


1 + X\ c 




(115) 

(116) 


Note that c is independent of Y. and depends only on 
the spacetime curvature at the light cone. (In Fig. ffl we 
see that all trajectories approach the light cone with the 
same slope c.) Integrating this we find that 


Q±(x,T\Y\) — G[x 0 (Y),Y] \x — x\ c \°, x~x\ c 

(117) 

We see that the decay of Q towards the light cone de¬ 
pends on c. When the metric is flat we have c = 1. This 
combines with a factor (x — Xi c ) arising elsewhere in the 
integrals to give a finite discontinuity at the light cone. 
However, if c < 1 the decay is faster and particles peel¬ 
ing off the light cone make no contribution at the light 
cone, giving a continuous stress-energy and particle cur¬ 
rent. Typically we find that c ~ 0.8 in our examples 
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when gravity comes into play. Therefore gravity makes 
the stress-energy tensor continuous at the light cone with¬ 
out the correlation between type I and type II particles 
that was required for this in flat spacetime. On the other 
hand, generically the stress-energy is not analytic at the 
light cone because of the non-integer power of \x — x\ c \, 
but we could have a situation where c -1 is an integer and 
then a suitable arrangement of the particle distributions 
could render an analytic stress-energy tensor. We have 
not been able to construct any particular example. 


F. Properties of the distribution function 


Vlasov matter is considered to be a good matter model 
both in Newtonian gravity and in general relativity be¬ 
cause it does not develop singularities in flat spacetime. 
For massive particles it has been proved that singulari¬ 
ties are also absent in the Vlasov-Poisson system even for 
large initial data] L4| , and for small data in the Einstein- 
Vlasov system fl5 ]. 

Even though it is clear that we cannot apply this sec¬ 
ond result to massless particles, it is interesting to com¬ 
pare the assumptions of that theorem with the properties 
of our solutions, because we can also construct solutions 
which are arbitrarily close to flat spacetime. The main 
assumptions of the theorem are that 1. the initial data 
has compact support in both momentum space and phys¬ 
ical space, 2. / is small, and in particular bounded, 3. / 
is C 1 and the metric is C 2 , and 4. the spacetime has a 
regular center. 

Using only type I particles, our metric solutions are 
generically only C 1 at the light cone, but we believe that 
correlating type I with type II particles solutions can be 
made C 2 there. We clearly have a regular center before 
the formation of the singularity, so we only need to an¬ 
alyze the boundedness of / and the compactness of its 
support at the initial time. 

Our CSS solutions are infinitely extended in space, but 
we can match initial data for a CSS solution for r < tq 
smoothly to initial data of compact support in space. As 
long as ro > x\ c (—t), the domain of dependence of the 
CSS part of the data includes the singularity. We then 
have a solution with compact support in space at the 
initial time which is nevertheless CSS in a central region 
including the singularity. 

We choose a function k(Y, Z) that has support only 
in a neighborhood of (lo,^o) not including Y = 0 or 
Z = 0. For given t < 0 and r, f(t,r,w,F ) then has 
support in momentum space only in a neighborhood of 
the two points ( w±(t , r, Yq, Zq), F±(t , r, Yq, Zq)). The two 
signs correspond to ingoing and outgoing particles, for 
x < x\ c . For x > Xic, only the positive sign applies. We 
find immediately that 


F±(t, r, V 0 , Z 0 ) = Z 0 t 2 Q ± 2 {x, Y 0 ), (118) 


which is finite for finite r. Because angular momentum 
is conserved along particle trajectories, we conclude that 


the support of / in angular momentum is always com¬ 
pact. 

w±(t, r, Yq, Zq) is the appropriate solution w of 


Y = 


x{x) 


1 + 


(rw) 




(119) 


for fixed r, t, Y = Y 0 and F = F±(t,r,Y 0 , Z 0 ). This is 
also finite. However, now w is not a constant of motion 
and therefore the support in w does not stay constant. 
The typical time evolution of the radial momentum w 
of a type I particle is the following: an ingoing particle 
starting with certain negative w accelerates towards the 
center, reaching a maximum modulus, and then decel¬ 
erates having a turning point and b eing ejecte d off the 
light cone with positive w. Eqs. (118) and (Ilf) to¬ 
gether show us that along lines of constant x, the radial 
momentum w of the particles that happen to be there 
at time t is independent of t. (The value of F of those 
particles even decays as F 2 . Note that at each t these 
are different particles.) Therefore the support in radial 
momentum remains compact during the evolution. 

We now show that in our CSS solutions f(t, r, w, F) can 
not be finite on the light cone r = x\ c {—t) at zero par¬ 
ticle momentum w = F = 0. Consider again a function 
k(Y. Z) that has support only in a small neighborhood 
of (Yo,Zo). Clearly if / in a CSS solution is bounded at 
one value of t < 0, it is bounded for all t < 0. Therefore 
we consider a fixed value t < 0. In order to find a limit 
in (r, w , F) space in which 


f{t,r,w,F) = ^k(Y,Z) 


( 120 ) 


blows up, we need to find a limit in which F —> 0 while 
Y —> Y 0 and Z —> Z 0 simultaneously. Consider therefore 
the limit 

F —* 0, w —> A(—f) _1 F 5 , x —> x\ c + BF 5 , 

_ (121) 
where c is the constant defined in (117), and A a nd B 
are real constants that will be determined. From (|llf 
we find that in this limit 


Y —» a(xi c ) (A 2 + x lc 2 ) 2 + A = Yq(A) 


( 122 ) 


and therefore 

Z ^ C[x 0 (Y 0 ),Y 0 ]- 2 \B\-i = Z 0 (A,B). (123) 


We can therefore always arrange the required limit for 
any values of lb >0 and Z 0 by a suitable choice of 
the constants A and B. We have therefore shown that 
/ is infinite on the light cone at zero momentum un¬ 
less it is identically zero. From the fact that the stress- 
energy tensor and the particle current are both finite if 
k(y) = f k(Y, Z)dZ exists, this blowup is not a physi¬ 
cal problem. Furthermore, for the ansatz of fc(F, Z) with 
compact support in Y and Z , bounded away from Y = 0 
and Z = 0, / is finite and can be made arbitrarily small 
(everywhere but at the light cone at zero momentum). 
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FIG. 10: Metric function a(x) and corresponding (dimensionless) mass function M(x). Note that the spacetime is not 
asymptotically flat, as we expected. The vertical lines give the position of the light cone x\ c = 1.3642. 
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FIG. 11: On the left we show the different iterations of the metric function a(x), starting with the flat case in the iteration 
i = 0. The convergence is fast and starting from i = 3 or i = 4 it is not possible to resolve different iterations in the figure. On 
the right, we show the decay of differences between successive iterations. We have defined Aid = cu — o»-i. The continuous 
line represents the 2-norm and the dotted line the oo-norm, both integrated between x = 0 and x = 2. 


V. CONCLUSIONS 

Type II critical phenomena, in which the black hole 
mass vanishes as a power of distance from the black hole 
threshold, have been found (in some region of parameter 
space) for almost all Einstein-matter systems in spherical 
symmetry. These include real and complex scalar fields 
with arbitrary potential terms, conformal couplings, and 
coupling to a Maxwell field, perfect fluids, sigma models 
and Yang-Mills fields. The only exception seems to be 
the spherically symmetric Einstein-Vlasov system. This 
raises the question what distinguishes this system from 
the other ones, and the wider question if the existence of 
type II critical phenomena is the rule or the exception. 

Critical phenomena at the black hole threshold require 
the existence of a critical solution. This is a solution that 
has precisely one unstable linear perturbation mode, with 
the additional property that a fully nonlinear evolution 
starting with a finite amplitude of this perturbation mode 
results in a black hole, while a finite amplitude of the 
opposite sign results in dispersion (or another outcome, 


such as a star or a naked singularity). In type II critical 
phenomena the critical solution is also self-similar, either 
continuously (CSS) or discretely (DSS). As a first step 
into understanding the absence of type II critical phe¬ 
nomena we have therefore investigated the existence and 
regularity of CSS solutions. 

The rest mass of the particles of the collisionless mat¬ 
ter introduces a scale into the coupled Vlasov and Ein¬ 
stein equations, which could be incompatible with ex¬ 
act self-similarity. However, there are many examples in 
other matter models where such a scale can be treated 
as a small perturbation in a class of solutions that are 
asymptotically self-similar. Type II critical phenomena 
are then not affected by the presence of the scale in the 
field equations. We have therefore focussed in this paper 
on spherically CSS solutions with a regular center of the 
Einstein-Vlasov system with massless particles, hoping 
to generalize these later to asymptotically CSS solutions 
with massive particles. 

The main result of this paper is the explicit construc¬ 
tion of a family of spherically symmetric CSS solutions 
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FIG. 12: Pressure function T r r (x). The 5-function case in 
flat space is represented with a dashed line; the Gaussian 
coupled to gravity is given by the continuous line and the 
Gaussian in flat spacetime is given by a dotted line which 
interpolates between the others. The box on the right shows 
an enlargement around the maximum. 



X 

FIG. 13: When Vlasov matter is coupled to gravity, the 
stress-energy tensor is continuous, but not C 1 . 


with massless particles that is parameterized by an ar¬ 
bitrary function of two conserved quantities k(Y, Z). By 
function counting we have constructed the most general 
spherically symmetric CSS solution with a regular center 
before the formation of the singularity, but there may be 
particular solutions that we have overlooked. We have 
also assumed the presence of certain cutoffs for very low 
and very large Y, Z , in order to avoid divergences. 

In order to make this result more transparent, we have 
also rederived the well-known general static spherically 
symmetric solution, which is parameterized by an arbi¬ 
trary function h(E, F) of the conserved particle energy E 
and angular momentum F. In both cases, static and CSS, 
we have first found the general solution for Vlasov test 
particles on a fixed spacetime of that symmetry, which 
is parameterized by a free function of three conserved 
quantities. We have then shown that when we demand 
that the resulting stress-energy tensor is compatible with 
the symmetry, the general solution depends only on two 
conserved quantities. In the static case, these are the ob¬ 


vious ones E and F, while the choice in the CSS case is 
much less obvious. 

Initial data for our solutions can be given compact 
support in momentum space, and can be truncated in 
space without affecting a central CSS region that includes 
the usual CSS singularity. However, the Vlasov function 
/(adiverges as —> 0 on the past light cone of 
the singularity. This divergence is integrable, so that 
both the stress-energy and the particle current are finite 
everywhere in spacetime except at the CSS singularity. 
Furthermore, our solutions can be constructed arbitrarily 
close to flat spacetime. 

The stress-energy tensor and particle current are less 
differentiable at the light cone than elsewhere. The 
stress-energy tensor of a generic test particle distribu¬ 
tion is discontinuous at the light cone, but by imposing 
a relation between k(Y, Z) and k(—Y,Z) it can be made 
C°. When the Vlasov matter is coupled to gravity, we 
gain one order of differentiability: the stress-energy ten¬ 
sor is now genetically C°, and can be made C 1 . The 
metric is genetically C 1 at the light cone, and can be 
made C 2 . (Naively, one would expect the metric to be 
two orders up from the stress-energy, but in polar-radial 
coordinates the two metric coefficients a and a can be de¬ 
termined from Einstein equations that contain only first 
derivatives of a and a.) 

Because our CSS solutions can be constructed with 
arbitrarily weak curvature, their curvature singularity is 
essentially kinematic: particles are “aimed” at the space- 
time point where the singularity will occur, rather than 
being focussed by gravity. In this context, the relation 
between regularity at the past light cone and the coupling 
to gravity is worth commenting on. Any spherically sym¬ 
metric CSS solution of a massless scalar test field that is 
regular at the center for t < 0 is necessarily singular at 
the past light cone. Coupling the scalar field to gravity, 
there exists an isolated strong field solution which is an¬ 
alytic both at the center and the light cone. The same 
is true for other field theories and a for a perfect fluid. 
By contrast, coupling collisionless matter to gravity adds 
one order of differentiability at the light cone to all mat¬ 
ter configurations, but no solution can be analytic at the 
light cone. This behavior is mathematically more similar 
to the behavior of the preferred scalar field solution at 
its future light cone. 

We shall consider the implications of our results for 
type II critical phenomena in detail elsewhere. How¬ 
ever, we have seen that the spacetime metric depends 
on the free function k(Y , Z) only through the integral 
k(Y) = J k{Y, Z)dZ. This means that there are infinitely 
many matter configurations that give rise to the same 
spacetime. It is clear that a similar result will hold for 
the linearized perturbations of these solutions. Therefore 
there will be an infinite number of linear perturbation 
modes with the same eigenvalue A. None of our solutions 
can therefore have a single growing mode. This seems to 
be related to the fact that Einstein-Vlasov is not a field 
theory. (While the phase space of the spherically sym- 
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metric scalar field consists of pairs (f>(r), 4>(r), the phase 
space of the spherically symmetric collisionless matter 
consists of functions /(r, w, F): this phase space is much 
bigger.) 

On the other hand, due to the fact that / is conserved 
along particles trajectories, we cannot expect to get a 
very close approach to a self-similar solution during the 
evolution of initial data with bounded /. We have seen 
that in the self-similar solutions / is unbounded for low 
momentum near the light cone, while numerical simula¬ 
tions typically work with finite /. 

Any, or both, of those two reasons could explain why 
the collapse simulations that have been carried out did 


not find any sign of type II critical phenomena. We must 
not forget, however, that those simulations worked with 
massive particles, while here we have assumed massless 
particles. 


Acknowledgments 

We would like to thank Hakan Andreasson and Alan 
Rendall for interesting conversations and comments. 
This research was funded in part by EPSRC grant 
GR/N10172. 


[1] M. W. Choptuik, Phys. Rev. Lett. 70 , 9 (1993). 

[2] C. Gundlach , Liv. Rev. Rel. 2 , 4 (1999), publ ished elec¬ 
tronically at pttp : //www.livingreviews.org . 

[3] G. Rein, A. D. Rendall, and J. Schaeffer, Phys. Rev. D 
58 , 044007 (1998). 

[4] I. Olabarrieta and M. W. Choptuik (2001), preprint gr- 
qc/0107076. 

[5] I. Olabarrieta, Master’s thesis, Univ. of British Columbia 
(2000), published as gr-qc/0012059. 

[6] A. Munier, J. R. Burgan, M. Feix, and E. Fijalkov, As- 
trophys. J. 236 , 970 (1980). 

[7] J. R. Gair, Class. Quant. Grav. (2001), to be published, 
preprint gr-qc/0110017. 

[8] J. Ehlers, in Relativity, Astrophysics and Cosmology , 


edited by W. Israel (Reidel, Dordrecht, 1973). 

[9] G. Rein, Habilitationsschrift, Univ. Munich (1995). 

[10] V. I. Arnold, Mathematical Methods of Classical Mechan¬ 
ics (Springer-Verlag, New York, 1989). 

[11] G. Rein, Math. Proc. Cambridge Philos. Soc. 115 , 559 
(1994). 

[12] J. Batt, W. Faltenbacher, and E. fforst, Arch. Rational 
Mech. Anal. 93 , 159 (1986). 

[13] J. Schaeffer, Comm. Math. Phys. 204, 313 (1999). 

[14] K. Pfaffelmoser, J. Differential Equations 95, 281 (1992). 

[15] G. Rein and A. D. Rendall, Comm. Math. Phys. 150 , 
561 (1992), erratum 176 , 475-478 (1996). 






